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ABSTRACT 

The  imaging  properties  of  a  uniaxial  anisotropic  slab  lens,  where 
the  dielectric  tensor  components  are  of  opposite  sign,  are  studied 
as  a  function  of  the  structure  parameters.  While  hypothetical 
parameters  yield  various  levels  of  performance,  a  design  principle 
to  achieve  good  sub-wavelength  resolution  is  suggested.  The 
anisotropy  can  be  implemented  with  a  metal-insulator  stack. 

The  influence  of  material  and  thickness  on  the  subwavelength  imaging 
performance  of  a  negative  dielectric  constant  slab  is 

studied.  Resonance  in  the  plane  wave  transfer  function  produces  a  high  spatial 
frequency  ripple  that  could  be  useful  in  fabricating  periodic 
structures.  A  cost  function  based  on  the  plane  wave  transfer  function 
provides  a  useful  metric  to  evaluate  the  planar  slab  lens  performance,  and 
using  this,  the  optimal  slab  dielectric  constant  can  be  determined. 

Prospects  for  a  lossless  negative  dielectric  constant  material  for 

optical  devices  are  studied.  Simulations  show  that  with  sufficient 

gain,  a  mixture  of  two  semiconductor  quantum  dots  can  produce  an  isotropic 

effective  dielectric  constant  that  is  lossless  and  negative.  Over 

length  scales  where  homogenization  is  meaningful,  this 

pennits  a  small-scale  optical  mode  volume  and  lossless  waveguides, 

major  goals  in  the  field  of  nanophotonics. 

We  show  quantitatively  the  influence  of  dipole  density  on  the 
resolving  capability  of  a  negative  refractive  index  slab.  A  numerical 
solution  for  scatter  with  a  model  that  represents  a  locally  periodic 
array  of  electric  and  magnetic  dipoles  provides  a  means  for  evaluating 
the  influence  of  periodicity.  A  Fourier  transfonn  of  the  spatial  field 
solution  gives  access  to  the  plane  wave  transfer  function  in  relation 
to  the  degree  of  granularity.  Throughout,  a 

quasi-static  homogenization  procedure  is  used  to  ensure  that  correctly 
matched  effective  medium  parameters  are  achieved.  The  results  provide 
a  measure  of  the  fabrication  technology  necessary  to  achieve 
substantially  subwavelength  resolution. 

We  suggest,  based  on  the  principle  of  causality,  and  for  a  material 
exhibiting  adjacent  absorptive  and  gain  resonances,  that  there  can  be 
an  intervening  frequency  where  perfect  imaging  is  in  concept 
possible.  At  this  frequency,  both  the  dielectric  constant  and  the 
permeability  are  negative,  leading  to  a  negative  refractive  index,  and 
there  is  no  loss.  In  such  a  material 

exhibiting  a  double  resonance,  the  gain  must  be  at  the  higher 
frequency.  Through  appropriate  tuning  of  the  refractive  index,  all 
propagating  and  evanescent  fields  from  the  object  could  then  in 
principle  be  reconstructed  at  the  image  plane,  subject  to  practical 
implementation  limits. 


The  work  reported  resulted  in  the  production  of  five  journal  papers. 
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Abstract 

The  imaging  properties  of  a  uniaxial  anisotropic  slab  lens,  where  the  dielectric  tensor  components 
are  of  opposite  sign,  are  studied  as  a  function  of  the  structure  parameters.  While  hypothetical 
parameters  yield  various  levels  of  performance,  a  design  principle  to  achieve  good  sub-wavelengtli 
resolution  is  suggested.  The  anisotropy  can  be  implemented  with  a  metal-insulator  stack. 

The  influence  of  material  and  thickness  on  the  subwavelength  imaging  performance  of  a  negative 
dielectric  constant  slab  is  studied.  Resonance  in  the  plane  wave  transfer  function  produces  a  high 
spatial  frequency  ripple  that  could  be  useful  in  fabricating  periodic  structures.  A  cost  function 
based  on  the  plane  wave  transfer  function  provides  a  useful  metric  to  evaluate  the  planar  slab  lens 
performance,  and  using  this,  the  optimal  slab  dielectric  constant  can  be  determined. 

Prospects  for  a  lossless  negative  dielectric  constant  material  for  optical  devices  are  studied.  Sim¬ 
ulations  show  that  with  sufficient  gain,  a  mixture  of  two  semiconductor  quantum  dots  can  produce 
an  isotropic  effective  dielectric  constant  that  is  lossless  and  negative.  Over  length  scales  where  ho¬ 
mogenization  is  meaningful,  this  permits  a  small-scale  optical  mode  volume  and  lossless  waveguides, 
major  goals  in  the  held  of  nanophotonics. 

We  show  quantitatively  the  influence  of  dipole  density  on  the  resolving  capability  of  a  negative 
refractive  index  slab.  A  numerical  solution  for  scatter  with  a  model  that  represents  a  locally  periodic 
array  of  electric  and  magnetic  dipoles  provides  a  means  for  evaluating  the  influence  of  periodicity. 

A  Fourier  transform  of  the  spatial  held  solution  gives  access  to  the  plane  wave  transfer  function  in 
relation  to  the  degree  of  granularity.  Throughout,  a  quasi-static  homogenization  procedure  is  used 
to  ensure  that  correctly  matched  effective  medium  parameters  are  achieved.  The  results  provide  a 
measure  of  the  fabrication  technology  necessary  to  achieve  substantially  subwavelength  resolution. 

We  suggest,  based  on  the  principle  of  causality,  and  for  a  material  exhibiting  adjacent  absorptive 
and  gain  resonances,  that  there  can  be  an  intervening  frequency  where  perfect  imaging  is  in  concept 
possible.  At  this  frequency,  both  the  dielectric  constant  and  the  permeability  are  negative,  leading 
to  a  negative  refractive  index,  and  there  is  no  loss.  In  such  a  material  exhibiting  a  double  resonance, 
the  gain  must  be  at  the  higher  frequency.  Through  appropriate  tuning  of  the  refractive  index,  all 
propagating  and  evanescent  fields  from  the  object  could  then  in  principle  be  reconstructed  at  the 
image  plane,  subject  to  practical  implementation  limits. 

The  work  reported  resulted  in  the  production  of  five  journal  papers. 

1  Summary  of  project 

Traditional  imaging  systems  or  lenses  operate  in  the  far  field  and  have  spatial  resolution  limited  by  the 
wavelength.  This  restriction  on  the  spot  size  occurs  because  of  the  diminishing  evanescent  fields  with 
increasing  distance  from  the  object,  and  the  image  is  formed  in  the  far-held.  This  means  that  imaging 
systems  where  resolution  is  important,  such  as  in  optical  lithography  for  semiconductor  devices,  must 
use  smaller  wavelengths  in  order  to  achieve  finer  features.  This  limitation  can  be  circumvented  by 
spatially  limiting  the  detection  area,  as  is  done  in  a  near-held  scanning  optical  microscope  (NSOM). 
While  the  NSOM  is  a  powerful  tool,  it  is  not  useful  for  many  applications  (lithography,  for  example), 
and  the  tip  must  be  mechanically  moved,  which  can  be  difficult  and  slow. 

If  it  were  possible  to  achieve  subwavelength  image  resolution  without  use  of  a  small  aperture  then 
the  impact  on  optics  would  be  enormous.  Imaging  elements  could  be  very  small,  the  resolution  of  optical 
microscopes  could  be  improved,  optical  interconnects  in  and  between  chips  could  be  practical,  optical 
storage  devices  could  contain  much  more  data,  optical  signal  processing  would  be  revolutionized,  and 
optical  lithography  could  be  used  to  achieve  much  small  chip  feature  sizes. 

Much  of  the  excitement  about  negative  index  materials  relates  to  two  possible  characteristics,  nega¬ 
tive  refraction,  which  allows  planar  lenses  and,  more  importantly,  with  small  loss,  growth  of  evanescent 
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fields,  which  can  facilitate  subwavelength  imaging.  As  there  are  no  known  materials  having  both  neg¬ 
ative  permittivity  and  permeability,  a  requirement  for  negative  index,  use  of  metamaterials  has  been 
explored.  The  difficulty  with  achieving  a  true  negative  index  material  that  provides  subwavelength 
resolution  is  the  need  for  both  high  dipole  density  and  low  loss,  and  the  need  for  a  magnetic  response  at 
the  wavelength.  As  we  have  shown,  any  loss  moves  operation  to  the  near  field  and  limits  the  resolution 
[1,  2]. 

Use  of  a  thin  metal  layer  “lens”  and  the  correct  polarization  can  produce  subwavelength  features  in 
the  near  field.  We  have  found  that  use  of  a  metal  film  stack  can  provide  significant  improvement  in  the 
evanescent  field  performance,  i.e.,  much  smaller  spot  size  [3].  We  have  done  a  preliminary  analysis  of 
a  metamaterial  film  structure  that  will  provide  the  best  resolution  of  any  imaging  device  that  does  not 
have  limited  spatial  support,  such  as  could  be  defined  by  a  small  aperture.  Our  plan  is  to  study  various 
structures  that  have  potential  for  subwavelength  imaging. 

2  Introduction 

Conventional  optics  suffers  from  the  diffraction  limit  because  only  the  propagating  components  of  the 
plane  wave  spectrum  are  available.  As  a  consequence,  the  spot  size  is  limited  to  about  one  wavelength. 
There  is  therefore  a  link  between  memory  density  and  wavelength,  for  example.  Information  about  the 
subwavelength  features  of  the  object  are  carried  by  the  evanescent  waves,  which  decay  exponentially. 
Subwavelength  imaging  is  of  fundamental  importance  in  fields  as  varied  as  semiconductor  device  lithog¬ 
raphy,  optical  memory  and  microscopy.  The  subwavelength  information  is  conveyed  in  the  evanescent 
(non-propagating)  portion  of  the  plane  wave  spectrum.  A  negative  refractive  index  slab  can  act  as  a  lens 
because  of  negative  refraction,  and  of  particular  interest,  can  amplify  the  evanescent  fields  and  produce 
images  with  subwavelength  features  [4,  5,  1,  2].  The  primary  features  which  negative  index  materials 
(NIMs)  offer  are:  negative  refraction,  evanescent  field  growth,  leading  to  the  prospect  of  subwavelength 
imaging,  and  dispersion.  These  characteristics,  individually  and  collectively,  are  important  for  imaging 
and  a  variety  of  other  optical  signal  processing  applications,  in  which  an  incident  field  (or  set  of  fields) 
is  transformed  in  some  manner  at  the  output  (such  as  to  a  subwavelength  spot).  While  building  nega¬ 
tive  index  materials  at  optical  wavelengths  will  be  very  challenging,  achieving  low-loss  structures  which 
perform  the  aforementioned  operations  is  possible  now. 

We  have  shown  that  a  layered  metal-insulator  stack  has  the  capability  of  subwavelength  imaging 
[6].  It  works  because  the  effective  dielectric  constant  is  anisotropic  and  hence  the  structure  converts 
evanescent  fields  from  the  source  (the  object)  to  propagating  wave  the  multi-layer  stack.  In  this  way, 
information  about  the  evanescent  fields  from  the  source  can  be  conveyed  to  the  image  plane. 


3  Results 


3.1  Sub- wavelength  imaging  opportunities  with  planar  uniaxial  anisotropic  lenses 


The  imaging  properties  of  a  uniaxial  anisotropic  slab  lens,  where  the  dielectric  tensor  components  are 
of  opposite  sign,  are  studied  as  a  function  of  the  structure  parameters.  While  hypothetical  parameters 
yield  various  levels  of  performance,  a  design  principle  to  achieve  good  sub-wavelength  resolution  is 
suggested.  The  anisotropy  can  be  implemented  with  a  metal-insulator  stack. 

Referring  to  Fig.  1,  we  assume  dielectric  constants  ex  =  e'x  +  je"  and  ez  =  e'z  +  jez,  with  time 
dependence  exp(ju;f).  For  the  TM  wave  (extraordinary  wave  with  Hy),  the  dispersion  relation  is  kxe~1  + 
k2e~l  =  kl  (and  in  free  space,  kx  +  k2z 0  =  fcg).  In  the  absence  of  loss,  if  ex  >  0  and  ez  <  0,  the  structure 
couples  the  incident  plane  wave  fields  (propagating  or  evanescent)  to  propagating  waves.  As  a  result, 
higher  spatial  frequency  Fourier  components  with  kx  >  ko  are  available  at  the  output  plane  and  hence 
for  subwavelength  imaging  [7,  8]. 

The  transmission  coefficient  for  Hy  is  given  by 


T(kx) 


1 

cos (kzd)  +  | (§0  +  §a)  sin (kzd) 


(1) 


where  Zs  =  kz/uexe o  and  Zq  =  kzo/Loeo  denote  the  characteristic  impedance  of  the  slab  and  the 
surrounding  medium,  respectively,  with  eo  the  free  space  permittivity. 
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Figure  1:  Anisotropic  slab  structure. 
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Figure  2:  Transmission  spectrum  for  (a)  lossy  or  optically  amplifying  media,  and  (b)  the  ez  — >  oo  lens.  The  parameters  are  A  =  700nm, 
d  =  400nm. 

For  the  deep  evanescent  spectrum,  the  magnitude  of  transmission  coefficient  remains  unchanged 
if  e"  and  e"  change  signs  simultaneously.  While  optical  gain  could  offset  the  adverse  impact  of  loss 
[9],  overcompensation  of  loss  will  detract  from  the  resolution.  Figure  2(a)  shows  an  example  where 
ex  =  1  —  j’0.05,  e2  =  —1  —  j0.05,  and  ex  =  1  +  j0.05,  ez  =  —1  +  j0.05,  have  essentially  identical 
\T(kx/ko)\  throughout  the  evanescent  spectrum  (kx/ko  >  1). 

Subject  to  low  loss,  tuning  exez  close  to  —1,  and  making  ex  less  positive  and  ez  more  negative  should 
improve  the  resolution,  and  some  loss  to  dampen  singularities  are  beneficial.  In  Fig.  3(a),  with  exjez 
fixed,  the  best  resolution  is  achieved,  as  anticipated,  with  ex  =  0.01  —  j'0.01  and  e2  =  —100,  for  which 
exez  is  closest  to  —1.  Figure  3(b)  shows  that,  with  exez  constant,  the  lens  with  the  lowest  value  of  \ex/ez\ 
(corresponding  to  ex  =  0.01  —  jO.Ol)  gives  the  best  resolution. 

With  ez  — >  oo  and  ex  — >  0,  a  pole  of  the  transfer  function  given  by 


occurs,  implying  that  the  pole  is  uniquely  determined  by  the  slab  thickness  d  at  a  fixed  frequency. 
Although  ez  — *  oo  and  ex  — >  0  do  not  make  a  perfect  lens  because  of  the  transmission  singularity,  these 
limits  together  with  small  slab  absorption  may  hold  promise  for  nearly  perfect  imaging.  The  underlying 
physical  mechanism  relies  on  the  fact  that  e2  — >  oo  reduces  sensitivity  to  the  transverse  wave  vector  kx , 
or  in  terms  of  the  “preferred  direction”  [10],  e2  — >  oo  eliminates  the  divergence  of  wave  components  by 
enforcing  the  “preferred  direction”  parallel  to  the  optical  axis,  and  ex  — ►  0  ensures  the  slab  is  electrically 
thin,  regardless  of  its  physical  thickness. 
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K  A 

d  =  n—  or  kzd  =  s/e^kod  =  mr,  n  =  1,2,  •••  (3) 

the  amplitude  of  the  transfer  function  becomes  unity,  and  under  this  resonant  condition,  the  lens  is 
perfect.  More  specifically,  kzd  =  2mr  results  in  a  “zero-phase  shift”  perfect  lens,  whereas  kzd  =  (2n+l)7r 
leads  to  a  ‘br-phase  shift”  perfect  lens.  Note  that  ez  — >  oo  fulfills  the  condition  of  the  canalization  regime 
[11,  12,  13].  A  small  deviation  of  ^/e^k^d  from  mr  introduces  a  pole,  thereby  compromising  the  system 
resolution.  Figure  2(b)  shows  the  transmission  coefficient  magnitude  for  a  perfect  lens  and  lenses 
operating  near  this  condition,  where  singularities  occur.  In  Fig.  3(c),  lenses  operating  in  the  multiple 
half- wavelength  regime  exhibit  similar  resolution  to  the  electrically  thin  lenses,  whereas  in  Fig.  3(d),  a 
small  deviation  from  the  perfect  condition,  say,  a  variation  of  thickness  d  leading  to  R e(kzd)  =  0.957T, 
significantly  impacts  the  field  profile  and  in  turn  the  resolution,  relative  to  the  electrically  thin  lens.  Our 
comparison  between  the  two  types  of  lenses,  as  given  in  Figs.  3(c)  and  3(d),  suggests  that  the  electrically 
thin  lens  may  be  superior  to  the  resonant  slab  lens  because  ex  — >  0  removes  the  strict  requirement  on 
the  physical  thickness.  The  resonant  perfect  lenses  in  Fig.  3(c)  do  not  restore  the  source  completely  in 
the  image  plane  because  of  the  mask  geometry,  but  the  field  on  the  left  interface  of  the  anisotropic  slab 
is  replicated  on  the  right  surface. 


X/A  X/A 

(a)  (b) 


X/A  X/A 


(c)  (d) 

Figure  3:  Double-slit  resolution  of  an  anisotropic  slab  lens,  based  on  finite  element  simulations,  with  various  material  parameters: 
(a)  ex/ez  is  constant;  (b)  exez  is  constant;  (c)  ex  — >  0,  ez  —>  oo  and  the  perfect  lens;  (d)  ex  — >  0,  ez  — >•  oo  and  deviation  from  the 
perfect  lens.  The  parameters  are:  A  =  700nm,  d  =  400nm. 


3.2  Imaging  performance  of  an  isotropic  negative  dielectric  constant  slab 

We  explore  the  near-field  imaging  properties  of  a  single  slab  of  material  with  isotropic  e  =  e'  +  je", 
with  e'  <  0.  This  condition  is  satisfied,  for  example,  by  metals  at  optical  frequencies.  It  is  also 
possible  to  synthesize  materials  that  have  these  properties,  such  as  mixtures  of  metals  and  dielectrics, 
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and  fabrication  can  be  relatively  straight  forward.  While  this  system  has  received  considerable  recent 
attention,  we  are  motivated  to  understand  what  resolution  range  is  accessible,  how  to  measure  resolution, 
and  potential  applications. 

We  consider  a  TM  held  and  a  free  space  background.  For  an  isotropic  slab  with  dielectric  constant 
e,  the  dispersion  relation  is  k2  +  k2  =  ek$,  where  in  free  space  =  w2/c2,  with  c  the  velocity.  Following 
the  exp  (jut)  time  convention,  the  transmission  coefficient  is 

(1  -  r2)exp{-jkzd ) 

1  —  r2exp(—j2kzd)  ’ 

where  r  is  the  reflection  coefficient  and  d  is  the  thickness  of  the  slab  [6].  Singularities  in  eq(4)  occur 
when  r2  =  exp  (j2kzd).  It  is  this  property  of  the  transmission  coefficient  approaching  infinity  that  we 
want  to  exploit  for  achieving  an  improved  image  resolution.  Also,  in  the  off-resonance  conditions,  we 
have  T(kx )  ~  (1  — r2)  exp(— kxd)  for  large  kx  and  real  e,  making  the  slab  thickness  a  dominant  parameter 
in  imaging  performance. 

A  perfect  lens  has  T(kx )  =  1  and,  generally,  the  resolution  can  be  expected  to  improve  as  T(kx ) 
increases  deep  in  the  evanescent  (free  space)  plane  wave  spectrum.  This  occurs  for  decreasing  slab 
thickness  d,  and  practical  considerations  drive  this  standoff  distance  between  the  object  and  image 
planes.  The  singularities,  which  dampen  with  loss,  are  undesirable  in  an  imaging  application. 

Within  the  object  spatial  support,  the  incident  field  is  defined  by  Hy  =  1,  i.e.,  the  incident  field  is 
defined  by  Hy(x)  =  p(x  —  a,w),  where  p(x  —  a,w)  is  a  unit  amplitude  pulse  centered  at  x  =  a  and  of 
width  w.  A  plane  wave  expansion  of  this  excitation  was  used  to  solve  for  Hy  in  each  region,  and  E 
was  determined  from  an  inverse  Fourier  transform  of  the  plane- wave  spectrum  for  Hy  multiplied  by  the 
impedance. 

Singularities  in  T(kx )  produce  dominant  image  plane  effects  corresponding  to  those  spatial  frequen¬ 
cies.  The  dominant  influence  of  the  singularity  causes  a  ripple  at  this  spatial  frequency.  Through  this 
mechanism,  by  control  of  the  complex  dielectric  constant,  it  is  in  concept  possible  to  generate  a  periodic 
spatial  ripple  with  arbitrarily  high  frequency.  This  could  be  useful  in  lithography  for  very  dense  (relative 
to  the  illuminating  wavelength)  periodic  systems  (photonic  crystals).  How  practical  this  is  rests  on  how 
easily  the  dielectric  constant  can  be  controlled,  in  particular,  the  imaginary  part. 

Passive  materials  will  have  some  loss,  and  we  explore  the  influence  of  e "  in  Fig.  4.  Fig.  4(a)  shows 
the  magnitude  of  the  transmission  spectrum  obtained  for  a  0.1  /mi  thick  slab  having  e1  =  —1.01.  Clearly, 
the  loss  suppresses  the  singularities.  With  an  object  Hy(x)  =  p(x,0.01A)  and  a  d  =  0.1  pm  thick  slab 
having  e'  =  —1.01  and  various  e”,  we  obtain  the  electric  field  images  of  Fig.  4(b).  The  full  width,  half 
maximum  (FWHM)  of  the  images  is:  0.265A  for  e  =  — 1.01— j'0.01;  0.235A  for  e  =  — 1.01  — j'0.02;  0.255A 
for  e  =  —1.01  —  j0.05;  and  0.295A  for  e  =  —1.01  —  j’0.1.  As  a  reference,  the  FWHM  after  propagation 
over  0.1  pm  in  free  space  is  found  to  be  0.415  A.  It  is  clear  from  this  example  that  the  best  performance 
is  not  obtained  when  the  loss  is  minimum.  This  leads  to  the  need  to  optimize  the  lens  parameters  for  the 
best  results.  A  well  chosen  cost  function  can  used  to  optimize  the  parameters  of  the  transfer  function. 


(4) 


e’  =  -1 .01 


e’ — 1.01 


(a) 


(b) 


Figure  4:  (a)The  magnitude  of  the  transmission  coefficient  for  a  0.1  /im  thick  slab  having  various  e  =  —1.01  —  je".  (b)  |E|.  The  line 
types  correspond  to:  e"  =  0.01  (solid)  e"  =  0.02  (dashed)  e"  =  0.05  (dash-dotted)  e"  =  0.1  (dotted)  The  free  space  wavelength  is 
A  =  700  nm. 
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where  F_1  is  the  inverse  Fourier  transform.  Equation  (5)  takes  into  consideration  both  the  magnitude 
as  well  as  the  phase  of  the  transmission  coefficient  and  has  been  used  for  our  results.  For  a  perfect 
lens,  (7  =  0.  With  A  =  700  nm,  we  searched  for  the  minimum  C(e,d )  for  a  series  of  d/X  using  a 
step  size  of  5d/X  ~  0.071.  We  found  the  optimum  e  at  each  d/ A,  based  on  the  minimum  of  (5),  and 
the  result  is  given  in  Fig.  5(a).  The  best  resolution  is  when  e'  is  slightly  less  than  —1.  Also,  e"  has 
been  optimized  for  minimum  cost.  Notice  that  the  optimal  e"  for  d/X  <  1  is  small,  but  non-zero. 
Fig.  5(b)  shows  the  FWHM  of  |E|  in  the  image  plane  for  the  optimal  e  shown  in  Fig.  5(a).  An  object  of 
width  0.01  A  was  used.  As  the  slab  thickness  increases,  the  lens  performance  deteriorates,  but  it  is  still 
substantially  better  than  the  free  space  case,  given  by  the  e  =  1  curve.  The  FWHM  is  subwavelength 
only  for  d/X  <  1,  making  this  the  useful  regime.  For  thicker  lenses,  a  larger  loss  {e")  gives  the  minimum 
cost  and  hence  smaller  FWHM.  This  loss  produces  a  propagating  component,  clearly  important  with 
increasing  thickness.  However,  this  large  loss  significantly  reduces  the  amplitude  of  the  electric  field  in 
the  image  plane.  Our  studies  thus  far  suggest  that  (5)  is  a  useful  metric. 


(a) 


(b) 


Figure  5:  (a)  Optimal  e  obtained  from  the  cost  function  (5),  for  various  slab  thickness  d.  (b)  Optimal  FWHM  for  |E|  as  a  function 
of  d.  The  d/X  scale  should  be  interpreted  as  various  d  with  A  =  700  nm.  The  object  in  (b)  has  width  0.01  A. 


3.3  Anisotropic  route  to  optical  circuits 

By  exploiting  the  optical  scattering  properties  of  a  single  slab  with  a  uniaxial  dielectric  tensor  having 
ek  =  ( k  =  x,y,z)  (as  in  Fig.  6),  we  find  that  for  the  transverse-magnetic  polarization  incident 

light  the  slab  can  play  a  role  of  (a)  a  distributed  transmission  line  (Fig.  6(a));  (b)  a  lumped  series 
inductor  (Fig.  6(b));  and  (c)  a  lumped  shunt  capacitor  (Fig.  6(c)).  In  contrast  to  the  realization  of 
circuit  nanoelements  in  the  optical  domain  using  nanoparticles  [14,  15],  the  present  approach  to  optical 
circuits  bears  two  compelling  virtues:  frequency-independence,  and  high  insensitivity  to  the  angle  of 
incidence  or  the  orientation  of  the  slab  relative  to  the  incident  wave.  Since  the  circuit  parameters  are 
uniquely  determined  by  the  dimensions  and  material  parameters  of  the  slab  (Fig.  6),  the  optical  circuit 
elements  can  effectively  operate  over  a  wide  range  of  frequencies  and  incidence  angles,  thereby  allowing 
us  to  export  and  transplant  various  microwave  and  RF  circuit  functional  designs  to  the  infra-red  and 
visible  frequency  regime  more  conveniently.  The  slab  is  assumed  infinitely  wide  in  the  derivation,  but 
practically  a  large  aspect  ratio  (width  to  length)  or  large  local  curvature  (relative  to  the  wavelength) 
suffices  the  circuit  approximation.  Finally,  we  note  that  the  desired  anisotropy  can  be  achieved  with  a 
metal-insulator  stack  or  wire  media  [16,  17,  8,  18]. 

An  optical  low-pass  filter  structure  is  depicted  in  Fig.  7(a).  We  used  both  Chebyshev  (equal  ripple) 
and  Butterworth  (maximally  flat)  polynomials  to  determine  the  L-C  values  having  a  normalized  cut-off 
frequency  ///o  =  1  [19].  Fig.  7(b)  shows  the  magnitude  of  the  transfer  function  for  Hy  for  both  normal 
( kx  =  0)  and  oblique  ( kx  =  0.8/co)  incidence.  The  frequency  roll-off  for  f  / fo  >  1  is  dictated  by  the  num¬ 
ber  of  sections.  Fig.  7(c)  shows  the  isotropic  layer  solution,  which  is  quite  similar  at  normal  incidence, 
but  deteriorates  rapidly  as  the  angle  of  incidence  is  moved  off  normal,  illustrating  the  robustness  to  the 
incident  angle  of  transmission  property  for  anisotropic  low-pass  filters  compared  to  the  isotropic  case. 

Impedance  matching/transforming  structures  provide  a  means  to  substantially  reduce  the  scattering 
cross-section.  Fig.  8(a)  illustrates  a  bilayer  L  —  C  section  that  can  eliminate  reflections  from  a  semi¬ 
infinite  Ag  layer.  Fig.  8(b)  shows  the  TM  ( Hy )  reflection  coefficient  magnitude  as  a  function  of  incidence 
angle  (kx/ko)  for  bilayer  matching  to  Ag.  The  two  types  of  capacitors  ( ez  — >  oo  and  ez  =  1)  are  shown, 
and  the  performance  with  the  e2  =  1  is  slightly  better.  For  comparison,  Fig.  8(b)  shows  an  isotropic 
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Figure  6:  (a)  Equivalent  transmission  line  in  the  optical  regime  from  envisioning  an  anisotropic  slab  when  ez  — ►  oo.  (b)  An  anisotropic 
slab  with  ez  — >  oo,  ex  — ►  0  could  play  the  role  of  an  optical  lumped  inductor  in  response  to  incident  TM  waves,  (c)  Equivalent  capacitor 
when  kzd  <C  1,  ex  1,  ez  — >  oo  (or  ez  =  1). 


(a) 


Cl  —  €()€xd±  L  2  —  /i0d2 


f/fo  f/fo 


Figure  7:  (a)  Synthesis  of  an  low-pass  filter  L-C  ladder  network,  (b)  Frequency  response  of  9-section  Chebyshev  (equal-ripple,  with 
film  thicknesses  li{i  =  1, . . . ,  9)  =  {28.1,  308.3,  37.2,  322.6,  37.6,  322.6,  37.2,  308.3,  28.1}(nm))  and  10-section  Butterworth  (maximally 
flat,  with  thicknesses  li(i  =  1, . . . ,  10)  =  {2.5,  360.8, 11.2,  708.0, 15.7,  784.9, 14.2,  561.9,  7.2, 124.3}  (nm))  low-pass  filters  with  normal 
and  oblique  incidence.  Each  anisotropic  film  assumes  ez  — >  oo,  and  ex  =  50  for  the  capacitor  and  ex  =  10— 3  for  the  inductor.  The 
wavelength  corresponding  to  fo  is  Ao  =  2500  nm.  (c)  Filters  with  the  parameters  in  (b)  and  isotropic  dielectric  constants:  e  =  50  for 
C  and  e  =  10-3  for  L. 


anti-reflection  coating,  and  this  shows  much  greater  sensitivity  to  angle  of  incidence.  Achieving  a 
low  reflection  for  a  wide  range  of  incidence  angles,  because  of  the  anisotropy  of  the  coating,  allows 
substantially  reduced  object  visibility.  To  illustrate  this,  Fig.  8(c)  shows  that  scattering  cross-section 
for  a  2-D  Ag  ellipse,  with  an  incident  TM  wave  along  the  semi-minor  axis.  Figs.  8(d)  and  (e)  show 
\Hy\  without  and  with  the  L  —  C  coating  (the  wave  is  incident  from  the  right).  The  result  for  a  bilayer 
having  e~  — >  oo  appears  superior.  The  improved  anti-reflection  properties  of  the  coating  could  also  find 


applications  in  efficient  solar  energy  collectors  and  other  electro-optical  devices  [20]. 
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Figure  8:  (a)  An  L  —  C  bilayer  anti-reflection  coating  on  Ag.  (b)  Bilayer  TM  reflection  coefficient.  Parameters:  A  =  400  nm, 
eAg  =  —3.77  +  z0.6747  [21],  ez  =  103,  ex  =  50.0  for  C,  ex  =  10— 3  for  L,  d\  =  124  nm,  d,2  =  3.0  nm.  The  result  for  an  isotropic 
quarter- wave  transformer  (e  =  3.0  and  length  92.6  nm),  on  a  spacer  layer  (e  =  45.4  and  length  14.9  nm)  provides  a  reference,  (c) 
Simulated  back-scattering  cross  section  of  the  scatter  with  and  without  the  bilayer  cover.  \Hy(x,  z)\  for  a  Ag  ellipse  under  TM  wave 
excitation  without  (d)  and  with  (e)  the  bilayer  cover. 


3.4  Lossless  negative  dielectric  constant  optical  material  based  on  a  semiconductor  quan¬ 
tum  dot  mixture 

We  show  that  a  metamaterial  with  a  negative  dielectric  constant  at  an  (optical)  frequency  where  the  loss 
goes  to  zero,  i.e. ,  e(u;o)  =  e' (uo)+ie" (u>o)  with  e"(u>o)  =  0,  can  be  devised  using  a  semiconductor  quantum 
dot  (QD)  mixture.  This  type  of  metamaterial  is  of  fundamental  importance  in  nanophotonics,  where  the 
down-scaling  of  high  density  optical  components  that  are  based  on  metal-insulator-metal  waveguiding, 
is  limited  by  the  high  losses  inside  the  metals  [22],  Exchanging  the  metals  by  a  semiconductor  QD 
mixture  with  a  similar  negative  dielectric  constant  but  no  losses,  would  allow  small  waveguide  mode 
volumes  and  hence  high  density  optical  elements  and  circuits  [23,  24]. 

Resonant  states  in  semiconductors  provide  a  source  for  negative  dielectric  constant,  provided  that  the 
dipole  strength  and  oscillator  density  is  adequate  to  offset  the  background.  Furthermore,  semiconductors 
offer  the  prospect  of  pumping,  either  optically  or  electrically,  to  achieve  a  gain  mechanism  that  may 
offset  the  loss.  Previously,  a  model  for  a  charged  QD  array  has  shown  a  negative  real  part  of  the 
dielectric  constant,  suggesting  isotropic  metal-like  properties  [25]. 

There  is  a  large  body  of  work  describing  the  scattering  of  light  from  QD  structures  that  provide 
electron-hole  (exciton)  confinement  in  three-dimensions  (for  example,  [26,  27,  28,  29,  30,  31,  32,  33, 
34]).  Models  typically  employ  an  effective  mass  approximation,  where  the  structure  is  assumed  large 
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relative  to  the  crystal  lattice  dimension.  In  the  case  of  a  spherical  QD  under  the  strong  confinement 
approximation,  we  assume  a  simplified  model  where  the  wave  function  can  be  written  as  the  product 
of  the  confined  electron  and  hole  functions.  This  makes  it  possible  to  write  an  analytic  solution  where 
the  energies  of  single  particles  can  be  summed  for  the  total  energy  [26,  28].  Assuming  also  that  the 
spherical  potential  well  has  infinite  barriers,  one  can  use  a  quantum  mechanical  derivation  to  arrive  at 
the  relation 

8e2  {fPaa  ~  l) 

eQD  €oc  Veo m*  lo1x  —  to2  —  i2ui^ 

where  is  the  matched  dielectric  constant  of  the  background  material  inside  the  QD,  e  is  the  electron 
charge,  eo  is  the  free  space  permittivity,  m*  is  the  effective  mass  of  the  electron,  uex  is  the  resonance 
angular  frequency  or  the  exciton,  7  is  the  phenomenological  damping  factor,  and  u  is  angular  frequency 

at  which  the  dielectric  constant  of  the  QD  is  calculated.  Also  in  (6),  pfg  is  the  steady  state  probability 
of  the  quantum  mechanical  exciton  system  being  in  the  ground  state.  Neglecting  thermally  excited 
electrons,  we  have  pfg  =  1  for  the  lossy  resonance  and  pfg  =  0  for  a  gain  resonance. 

The  dielectric  constant  e(u>)  for  the  ensemble  of  the  QDs  can  be  described  by  the  dielectric  theory 
of  Maxwell- Garnett  (MG)  [35,  36,  37,  38],  which  for  two  quantum  dot  species  can  be  written  as 


(6) 


e((j)  -  €i(u)  £QDlM  -  €i((j)  CQD2M  “  ti(u) 

-  ——  X\ - - -  X“2 - - - 

e(ui)  +  2ei(iv)  €qdi(w)  +  2  ti(u)  €QD2(w)  +  2ei{ui) 


(7) 


where  ay  is  the  volume  fraction  of  the  i-th  QD  species  and  et  is  the  background  dielectric  constant.  It 
is  assumed  that  the  QDs  are  spherical  and  that  macroscopic  dielectric  constants  can  be  applied.  Equa¬ 
tion  (7)  provides  an  approximate  means  to  account  for  polarization  charge  from  a  mismatch  between 
the  QD  and  background  dielectric  constant.  Solving  (7)  for  e  gives 


_  £i(  1  +  2 Q) 

~  1  -Q  ' 


We  focus  on  CdSe  QDs  and  use  Equation  (6)  to  find  the  QD  dielectric  constant,  and  then  (7)  and 
(8)  to  determine  the  homogenized  dielectric  constant  for  the  QD  mixture  in  a  background.  For  CdSe, 
m*e/nriQ  =  0.1,  where  mo  is  the  electron  rest  mass,  and  =  10.2.  Also,  the  exciton  resonance  energy 
can  be  calculated  for  CdSe  QDs  with  different  radii.  For  example,  ftwex  =  3.23  eV  for  R  =  1.5  nm, 
ftwex  =  2.55  eV  for  R  =  2  nm,  and  fhvex  =  1.91  eV  for  R  =  4  nm.  Based  on  past  fabrication  techniques 
that  produced  CdSe  QDs  with  linewidths  of  fry  =  0.04  eV  at  room  temperature  and  fry  =  0.02  eV  at 
10  K,  (and  were  limited  by  the  variation  in  dimension  of  the  QDs)  [33],  we  optimistically  assume  that 
state  of  the  art  QD  fabrication  technique  can  reduce  the  linewidths  to  around  fry  =  0.01  eV,  or  lower. 

Panel  (a)  of  Figure  9  gives  €qd  for  the  CdSe  QD  for  several  fry  values.  Notice  that  fry  <  0.02  is 
necessary  to  achieve  a  negative  dielectric  constant.  The  dependency  of  €qd  on  R  is  shown  in  panel  (b) 
of  Figure  9.  From  this  figure,  it  is  evident  that  QDs  with  higher  exciton  energies  and  smaller  R  values 
offer  more  negative  dielectric  constants. 

We  require  e"  =  0  at  some  frequency  in  order  to  achieve  lossless  surface  plasmons  and  hence  lossless 
waveguide  modes.  To  investigate  prospects  for  obtaining  isotropic  material,  consider  distributions  of 
two  types  of  QDs,  one  providing  an  absorptive  resonance  and  the  other  gain.  The  gain  resonance  is 
assumed  to  be  at  ftwex  +  hS.  where  hwex  is  the  loss  resonance  and  hS  the  difference  energy  achieved 
through  a  small  reduction  in  dot  size.  Equal  lossy  and  gain  QD  fill  fractions  are  assumed,  i.e. ,  x\  =  xy. 
Panel  (a)  of  Figure  10  gives  example  results  that  indicate  that  a  condition  of  e"  =  0  and  e'  <  —  1  can  be 
achieved  with  moderate  to  high  fill  fractions  and  a  sufficiently  low  linewidth  fry  =  0.007  eV.  Panel  (b) 
of  Figure  10  shows  a  comparison  between  each  of  the  plots  of  the  effective  dielectric  constant  in  panel 
(a),  obtained  using  the  MG  mixing  rule,  and  the  effective  dielectric  constant  obtained  via  a  different 
effective  medium  parameter  extraction  procedure.  In  this  procedure,  the  effective  medium  parameters 
are  extracted  from  a  full- wave  electromagnetic  simulation  of  the  transmission  through  a  slab  comprising 
periodic  inclusions  of  QDs  [39].  While  for  a  low  filling  factor  (x\  =  X2  =  0.1),  the  effective  dielectric 
constant  values  are  very  close  to  those  obtained  by  the  MG  mixing  rule,  the  effective  dielectric  constant 
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(a) 


(b) 


Figure  9:  (a)  Dielectric  constant  for  a  CdSe  QD  with  R  =  2  nm  and  hy  =  0.01  eV  (squares),  0.02  eV  (no  symbol),  and  0.04  eV 
(circles),  (b)  Dielectric  constant  for  a  CdSe  QD  with  varying  R  and  foy  =  0.02  eV:  R  =  1.5  nm  (no  symbol),  R  =  2  nm  (circles), 
R  =  4  nm  (squares). 


values  for  a  higher  filling  factor  (aq  =  X2  =  0.2)  differ,  to  some  extent,  from  the  corresponding  MG 
mixing  rule  values.  We  conclude  that  the  MG  mixing  rule  can  give  a  good  qualitative  approximation 
to  behavior  of  the  QD  mixture,  which  helped  in  the  initial  prediction  that  the  condition  of  e"  =  0  and 
e'  <  —  1  could  be  met  for  a  loss-gain  QD  mixture.  However  for  future  experimental  pursuits,  where 
accurate  prediction  of  the  values  of  e'  is  important,  a  numerical  effective  medium  parameters  extraction 
procedure,  such  as  the  one  we  compared  the  MG  results  with  in  this  work,  will  be  required. 


Figure  10:  (a)  Effective  dielectric  constant  for  R  «  2  nm  CdSe  QDs  with  =  10  and  HS  =  0.06  eV,  with  x\  =  X2  =  0.2  (no  symbols) 
and  xi  =  X2  =  0.1  (circles),  assuming  hy  =  0.007  eV.  The  gain  is  at  the  higher  frequency,  (b)  Comparison  of  each  of  the  plots  in  (a) 
that  were  obtained  using  Maxwell- Garnett  (M-G)  mixing  rule,  with  plot  obtained  via  a  full-wave  (F-W)  effective  medium  parameter 
extraction  procedure. 


3.5  The  influence  of  granularity  in  a  negative  refractive  index  lens 

We  demonstrate  the  influence  of  the  lattice  geometry  or  the  degree  of  granularity  on  the  image  resolution 
of  a  negative  index  lens.  This  provides  a  quantitative  measure  of  the  necessary  dipole  density  for 
achieving  substantially  subwavelength  resolution.  Our  results  thus  place  experimental  programs  to 
achieve  enhanced  resolution  at  optical  frequencies  in  perspective. 

A  simple  measure  of  the  influence  of  granularity  is  the  finite  difference  solution  of  Maxwell’s  equa¬ 
tions,  where  a  point-wise  representation  is  used.  A  numerical  dispersion  relation  can  be  developed  for 
the  discrete  problem  as  [40] 


smc 


kx  A X  , 

*mi  + 


(  kz  Az 
Vh  A0 


=  MG 


(9) 


where  Ax  and  Az  are  the  mesh  dimensions  in  the  x-  and  ^-directions,  respectively,  sinc(x)  =  (sinx)/x, 
/r  is  the  relative  permeability,  and  e  is  the  dielectric.  The  exact  analytical  dispersion  relation  is  given 
by  kxa  +  =  k^fie.  We  set  kx  =  kxa  and  use  (9)  to  find  the  kz  =  kza  +  5kz  for  the  discrete  problem, 

where  5kz  is  the  error  due  to  the  discrete  representation.  Figure  11  shows  the  normalized  discretization 
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(a)  (b) 


Figure  11:  Numerical  dispersion  error  for  a  square  finite  difference  mesh  of  dimension  A.  (a)  \5(kz/ko)\  as  a  function  of  kx/ko  for: 
A/Ao  =  1/10  (dashed  line);  A/Ao  =  1/20  (dash-dotted  line);  A/Ao  =  1/30  (solid  line);  A/Ao  =  1/50  (dotted  line),  (b)  \S(kz/ko)\  as 
a  function  of  A/Aq  for:  kx/ko  =  0  (dashed  line);  kx/ko  =  2  (dash-dotted  line);  kx/ko  =  4  (solid  line);  kx/ko  =  6  (dotted  line). 


Image  Plane 


d/2 


Figure  12:  Discrete  negative  refractive  index  slab,  where  the  electric  and  magnetic  dipoles  are  represented  by  a  periodic  structure 
with  material  having  effective  e  =  —  1  and  fi  =  —  1. 


error  \5(kz/ko)\  as  a  function  of  kx/k o  and  A/Ao-  The  fields  are  evanescent  for  kx/ko  >  1.  Note  that 
the  error  for  the  evanescent  fields  grows  dramatically  with  kx/k o  and  A/Ao-  A  common  guideline  for 
computational  domains  of  a  few  wavelengths  is  that  at  least  20  grid  points  per  wavelength  are  needed 
to  capture  the  propagating  held  phase  front.  Taking  A/Ao  =  0.05  in  Fig.  11(b)  for  the  kx/ko  =  0 
curve,  a  comparable  error  for  the  kx/ko  =  6  plane  wave  requires  A/Ao  ~  2  x  10-3  (500  grid  points  per 
wavelength) . 

To  understand  the  impact  of  granularity  in  a  negative  refractive  index  slab,  we  propose  a  model  that 
makes  for  a  tractable  numerical  problem.  Figure  12  shows  the  model  we  use  for  a  discrete  metamaterial 
slab,  where  circular  cylinders  on  a  square  lattice  having  parameters  £2  and  /j,2,  selected  appropriately, 
produce  homogenized  parameters  e  and  fj,.  Each  cylinder  in  Fig.  12  is  a  model  for  coincident  electric  and 
magnetic  dipoles.  We  do  not  consider  how  these  parameters  may  be  achieved  physically  (for  example, 
through  some  arrangement  of  conductors),  only  that  some  £2  and  /J2  give  rise  to  e  =  —1  and  fi  =  —  1.  The 
incident  held  is  excited  in  a  very  small  region  in  the  object  plane  above  the  slab.  Varying  the  density  of 
the  cylinders  allows  us  to  quantitatively  appraise  the  impact  of  granularity  on  image  resolution.  In  our 
calculations  we  assume  a  TM-polarized  held  (H  parallel  to  the  cylinders).  Assigning  effective  medium 
parameters  requires  operation  near  the  center  of  the  Brillouin  zone,  i.e.,  with  sufficient  dipole  density. 
We  apply  considerable  care  in  forming  homogenized  parameters  that  result  in  a  good  approximation  to 
n  =  —  1,  because  our  study  requires  the  influence  of  material  mismatch  to  be  small  relative  to  that  for 
granularity.  We  should  note  that  the  periodic  system  of  Fig.  12  has  discrete  Bloch  modes  (in  kx)  at 
each  A’oy/ue  [41],  and  that  we  solve  for  the  scattered  helds  without  enforcing  periodicity  to  capture  the 
imaging  performance. 

For  a  TM-polarized  incident  wave,  we  determine  the  effective  //  using  /j  =  +  /xi(l  —  /),  where 

/  is  the  hll  fraction,  i.e.,  the  ratio  of  the  circle  to  the  unit  cell  areas  in  Fig.  12.  The  continuity  of  the 
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(a)  (b) 

Figure  13:  (a)  Transfer  function  for  a  homogeneous  slab  lens:  d  =  0.5Ao  and  e"  /e*  =  I  O'  3  (solid  line);  d  =  0.5Ao  and  e"/e'  =  1 0 ~  2 
(dashed  line);  d  =  Ao  and  e" /e'  =  10-3  (dotted  line).  The  symbols  are  the  results  using  the  plane  wave  transfer  function  for  each 
case,  (b)  The  transfer  function  for  the  discrete  slab  lens  (Fig.  12)  for  d  =  0.5Ao,  a  fill  fraction  of  0.5,  ei  =  1,  £2  =  —1.377,  // 1  =  I 
and  //2  =  —3:  30  cells  per  Aq  (dotted  line);  40  cells  per  Aq  (dashed  line);  50  cells  per  Aq  (solid  line). 


magnetic  field  at  the  surface  of  the  cylinder  makes  this  a  satisfactory  choice.  With  fill  fraction  /  =  0.5, 
/12  =  —3  and  ji\  =  1  result  in  \i  =  —1,  and  we  assume  this  case.  To  calculate  the  effective  e,  we  use 
a  numerical  electrostatic  model  [42].  We  solve  Laplace’s  equation  for  the  electrostatic  problem  of  a 
cylindrical  rod  in  a  square  cell,  as  in  Fig.  12,  under  the  assumption  that  the  cell  size  is  small  relative  to 
the  optical  wavelength.  Considering  the  case  with  fill  fraction  0.5,  we  choose  €2  =  —1-377  and  e\  =  1  to 
produce  e  =  —  1 . 

A  finite  element  method  was  used  to  solve  the  TM  wave  equation 

V  •  ( e~1’VHy(x ,  z ))  +  kQfiHy(x,  z)  =  0  (10) 

for  the  two-dimensional  scattered  magnetic  field  Hy(x,z),  where  ko  =  lx/c  is  the  wave  number  (cu  is 
the  circular  frequency  and  c  is  the  speed  of  light).  The  computational  domain  was  terminated  in  a 
Sonmierfeld-type  radiation  boundary  condition.  The  slab  width  will  affect  the  image  resolution  [43],  as 
will  truncation  of  the  data  set  in  space.  We  made  the  computation  domain  as  large  as  possible  and 
established  an  adequate  mesh  density,  given  the  constraint  on  the  number  of  unknowns  presented  by 
computer  memory.  In  the  simulation  of  the  discrete  slab,  a  sufficiently  large  discrete  region  (of  thickness 
d,  referring  to  Fig.  12)  was  placed  between  regions  having  continuous  material.  The  reason  for  this  is 
the  granular  region  requires  a  much  high  finite  element  mesh  density,  and  this  allows  a  compromise  to 
extract  results  of  sufficient  accuracy  to  allow  conclusions  on  the  impact  of  granularity.  In  addition,  to 
alleviate  the  influence  of  a  finite  slab,  we  placed  a  radiation  boundary,  assuming  a  homogeneous  slab, 
along  the  regions  of  the  computational  boundary  defining  the  ends  of  the  slab,  making  the  slab  appear 
to  be  (approximately)  infinite  in  extent.  The  field  in  the  image  plane  ( Hyi{x ))  and  the  object  plane 
(Hyo(  x))  were  then  Fourier  transformed  and  the  ratio  formed  to  determine  the  complex  transmission 
coefficient  as  a  function  of  the  transverse  spatial  spectral  variable, 


r(kx) 


Hyi{kx) 

Hyoikx) 


(11) 


where  kz  =  ^/k^  —  k%.  Evaluating  the  performance  of  r(kx)  for  the  evanescent  fields  (kx  >  ko)  provides 
an  important  measure  on  possible  subwavelength  imaging  performance. 

To  establish  the  accuracy  of  our  numerical  approach,  we  compare  the  finite  element  solution  for 
a  homogeneous  slab  lens  with  the  exact  plane  wave  solution  [2].  Figure  13(a)  shows  various  transfer 
functions  for  a  homogeneous  slab  lens  with  assumed  continuous  material  parameters  (£2  =  —  e'  +  ie", 
with  e'  =  1,  and  /i2  =  —1)  and  a  free  space  background  (ei  =  1  and  =  1),  calculated  using  the  plane 
wave  transfer  function  [2]  and  by  the  finite  element  method,  with  a  spatial  Fourier  transform  of  the 
fields  in  the  object  and  image  planes.  Notice  that  we  are  able  to  capture  the  low-pass  characteristics 
with  the  finite  element  solution,  and  that  these  vary,  depending  on  the  degree  of  loss.  This  data  provides 
a  reference  for  the  potential  accuracy  in  solving  the  discrete  structure  problem  of  Fig.  12.  The  finite 
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Figure  14:  The  resolution  R\r\  as  a  function  of  A/Ao  for  d  =  0.5Ao,  a  fill  fraction  of  0.5,  e±  =  1,  €2  =  —1.377,  fi i  =  1,  and  fi2  =  —3. 

element  results  for  a  discrete  slab  having  30,  40,  and  50  cells  per  free  space  wavelength  (Ao)  are  given 
in  Fig.  13(b),  where  effective  e  =  —1  and  /i  =  —1,  as  determined  using  the  homogenization  procedure 
we  have  described.  Increasing  granularity  in  this  range  reduces  the  plane  wave  spectrum  bandwidth, 
as  does  dissipative  loss.  It  appears  that  the  required  dipole  density  can  be  viewed  as  a  density  relative 
to  the  wavelength  in  the  transverse  direction,  \x  =  2ir/kx.  In  order  to  maintain  an  adequate  finite 
element  mesh  density,  for  an  accurate  solution,  we  were  limited  to  a  cylinder  cell  density  of  about  50 
cells  per  Ao-  Subwavelength  resolution  is  related  to  the  bandwidth  of  the  plane  wave  spectrum,  and  it  is 
clear  from  Fig.  13(b)  that  a  very  high  density  is  needed  to  correctly  capture  a  significant  portion  of  the 
evanescent  field  spectrum.  From  Fig.  13(b),  the  |r(fcx//co)|/v/2  point  with  50  cells  per  Ao  corresponds 
to  kx/k0  2.  The  finite  difference  dispersion  error  results  of  Fig.  11  indicate  that  the  error  in  kz ,  5kz, 
at  this  density  and  for  kx/ko  =  2,  is  similar  to  that  with  a  density  of  20  mesh  points  per  wavelength  for 
the  propagating  waves.  It  thus  appears  that  the  finite  difference  numerical  dispersion  error  model  of  (9) 
may  indeed  be  a  useful  measure  to  determine  the  necessary  level  of  granularity  in  a  negative  index  slab. 
We  define  the  resolution  i?|Ti  =  kxm/ko ,  where  kxm/ko  is  the  value  at  |r|  =  l/\/2-  Fig-  14  shows  our 
calculated  results  for  the  resolution  as  a  function  of  A/Ao-  Appreciable  evanescent  field  content  clearly 
requires  a  very  high  structure  density  ((A/Ao)-1). 

The  Poynting  vector  results  for  the  cases  of  Fig.  13  provide  image  information.  Figure  15  shows 
Sz ,  the  z-component  of  the  Poynting  vector,  for  the  homogeneous  slab  (Fig.  15(a))  and  the  discrete 
slab  (Fig.  15(b)),  normalized  by  the  maximum  of  the  ^-component  of  the  Poynting  vector  in  the  object 
plane.  The  image  resolution  improves  with  reducing  loss  and  increasing  dipole  (cylinder)  density.  This 
finding  is  consistent  to  the  results  of  the  transfer  functions  shown  in  Fig.  13.  Notice  that  the  Poynting 
vector  in  the  image  plane  shows  negative  power  density  regions.  These  are  associated  with  vortices  that 
we  found  in  negative  index  slabs  with  small  amount  of  loss  [2],  and  they  occur  here  because  of  the  level 
of  granularity. 

The  quantitative  numerical  studies  of  the  influence  of  discreteness  in  a  negative  index  slab  lens 
indicate  that  a  very  high  density  of  electric  and  magnetic  dipoles  will  be  necessary  to  achieve  substan¬ 
tially  subwavelength  resolution.  This  may  preclude  the  use  of  lithography  as  a  means  to  define  metallic 
dipoles  at  optical  wavelengths,  should  subwavelength  feature  resolution  be  the  goal,  and  suggests  that 
alternate  strategies  be  considered. 

3.6  A  perfect  lens  material  condition  from  adjacent  absorptive  and  gain  resonances 

Perhaps  the  most  important  property  of  a  slab  composed  of  a  negative  refractive  index  (or  left-handed) 
material  [4]  is  its  ability  to  amplify  evanescent  plane  waves  [5],  subject  to  loss,  mismatch  in  the  real 
part  of  the  refractive  index,  and  the  dipole  density  [44,  1,  45].  A  sufficient  requirement  for  a  negative 
refractive  index  is  that  the  dielectric  constant  and  the  permeability  be  simultaneously  negative.  While 
naturally  occurring  negative  index  material  is  not  available,  conducting  electric  and  magnetic  dipoles 
(metamaterials)  operating  in  the  effective  medium  limit  have  been  shown  to  exhibit  a  negative  refractive 
index  [46].  Of  note,  zero  loss  at  the  operating  frequency  would  potentially  allow  perfect  imaging,  the 
prospect  of  which  is  of  enormous  importance  in  microscopy,  semiconductor  device  processing,  and  optical 
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Figure  15:  The  z-component  of  the  Poynting  vector  in  the  image  plane  normalized  by  the  maximum  of  ^-component  of  the  Poynting 
vector  for  the  incident  field,  (a)  Homogeneous  slab  (Fig.  13(a)):  d  =  0.5Ao  and  e" /e'  =  10— 3  (solid  line);  d  =  0.5Ao  and  e" /e'  =  10-^ 
(dashed  line);  d  =  Ao  and  e" /e'  =  10“ 3  (dotted  line),  (b)  Discrete  slab  (Fig.  13(b))  with  d  =  0.5Ao,  a  fill  fraction  of  0.5,  ei  =  1, 
62  =  —1.377,  ^i  =  l  and  fi2  —  —  3:  30  cells  per  Ao  (dotted  line);  40  cells  per  Ao  (dashed  line);  50  cells  per  Ao  (solid  line).  The 
dash-dotted  lines  in  (a)  and  (b)  show  the  ^-component  of  the  incident  Poynting  vector  in  the  object  plane. 


memory.  There  has  been  interest  in  using  gain  to  offset  loss,  which  has  led  to  the  fundamental  question 
as  to  whether  it  is  possible  to  achieve  a  zero  loss  operating  condition  with  the  dielectric  constant  and 
the  permeability  both  negative  [47].  We  show  for  a  causal  system  described  by  multiple  resonances  that 
adjacent  gain  and  loss  resonances  provide  the  possibility  of  achieving  the  material  property  characteristic 
of  a  perfect  lens  [5],  and  a  resolution  subject  to  the  geometrical  arrangement. 

The  Kramers-Kronig  relations,  which  describe  the  spectral  relationship  between  the  real  and  imag¬ 
inary  parts  of,  for  example,  the  dielectric  constant,  are  widely  used  [48].  These  relations  form  a  Hilbert 
transform  pair  that  results  from  the  application  of  the  residue  theorem  of  Cauchy  [49]  in  the  complex 
u  =  u>'  +  iu"  plane,  and  are  given  by 
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where  (  =  +  iQ"  and  (  =  {//,  e,n2},  with  /i  the  permeability,  e  the  dielectric  constant,  and  n2  the 

square  of  the  refractive  index,  and  the  integral  is  interpreted  as  a  principal  value,  with  the  pole  residue 
evaluated.  The  requirement  for  these  relations  to  hold  is  that,  for  exp(— iiot),  £( to )  be  analytic  in 
the  upper  half  ui  plane.  Satisfaction  of  the  Kramers-Kronig  relations  guarantees  causality,  as  does  an 
analytic  e(u)  in  the  upper  half  plane  in  complex  lo  space.  An  assumed  solution  to  (12)  leads  to  our 
claim  that  a  perfect  lens  is  in  principle  possible. 

Equation  (12)  is  based  on  the  (electric  and  magnetic)  response  being  zero  before  some  retarded, 
causal  time,  tr  =  l/c,  where  l  is  the  distance  between  the  source  and  dipole  position  and  c  is  the  speed 
of  light  in  vacuum,  regardless  of  the  material  parameters  [50].  Use  of  group  velocity  in  establishing 
the  causal  time,  i.e. ,  energy  transport,  is  imprecise  (because  of  the  wave  packet  basis  in  the  Fourier 
domain),  and  it  is  erroneous  to  use  the  phase  velocity,  which  is  opposite  in  direction  to  the  Poynting 
vector  (the  direction  of  power  flow)  for  a  negative  index  material.  The  Kramers-Kronig  relations  do 
not  convey  the  causal  time,  just  that  there  is  a  local  time  ( t  =  0)  prior  to  which  the  dipole  moment  is 
zero.  They  can  thus  be  applied  to  any  realizable  system,  with  or  without  gain  (stimulated  emission), 
and  with  positive  or  negative  refractive  index,  but  not  one  with  upper  half-plane  poles. 

From  a  classical  oscillator  development  [48]  or  a  quantum  mechanical  derivation  [51],  the  electric 
susceptibility  can  be  written  in  the  form 
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where  e  it  the  electron  charge,  eo  is  the  free  space  permittivity,  m  is  the  mass,  A Nj  is  the  oscillator 
density  for  the  j-th  resonance  (A  Nj  >  0  corresponds  to  loss  and  A  Nj  <  0  to  gain),  fj  is  a  normalized 


15 


oscillator  strength  QW  fj  =  1),  and  7 j  is  the  linewidth  (7-  1  is  the  photon  lifetime).  The  linewidth  is 
broadened  by  inelastic  collisions  (with  phonons)  and  elastic  phase-destroying  collisions.  We  consider 
that  the  electric  and  magnetic  response  (for  a  metamaterial  operating  in  the  effective  medium  limit 
where  homogenization  is  appropriate)  can  be  represented  by  responses  of  the  form  of  (13),  as  can  a  gain 
mechanism.  The  poles  of  (13)  occur  at 

cu  =  ±^(4^-7 (14) 

where  7 >  0  for  both  loss  (absorption)  and  gain  (stimulated  emission).  Any  x(cu)  that  can  be  expressed 
in  the  form  of  (13)  is  both  causal  and  analytic  in  the  upper  half  of  the  complex  to  plane. 

Consider  now  the  case  of  two  resonances  in  (13),  one  with  loss  and  one  with  gain.  The  fundamental 
question  of  interest  is  whether  it  is  possible  to  have  x"{ljJp)  =  0  (e”(Lop)  =  0  and  ^”(cop)  =  0)  at  a 
frequency  where  x'  <  —  1  (both  e'(ujp)  <  0  and  //( ojp )  <  0).  We  write  the  (electric  or  magnetic) 
susceptibility  as 


°i  .  «2  nr, 

x  =  — - 2 — ; —  +  — - 2 — : — '  (15) 

uf  —  iv~  —  iuo  71  uj  2  —  co  —  1^X2 

As  a  simplification,  let  71  =  72  =  7,  and  a\  =  —(12  =  a,  where  the  latter  condition  forces  one  of  the 
resonances  to  offer  gain.  Using  these  parameters  and  solving  (15)  for  =  Xi  +  X2  =  0  leads  to 


LU  =  CU p 


Letting  27s  =  (uq  —  uq)2  gives 
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In  (17),  x'i^p)  <  0  when  a  >  0  (loss  at  uq  and  gain  at  u q),  assuming  u  1  <  uq. 

The  (electric  or  magnetic)  susceptibility  in  (15)  is  causal  and  hence  is  a  physically  meaningful 
description  of  a  metamaterial.  As  x'{up)  can  be  negative  with  the  appropriate  choice  of  parameters, 
it  is  possible  to  have  both  7/  and  e'  negative  at  a  frequency  where  both  n"  and  e"  are  zero.  Figure  16 
shows  the  case  of  (15),  with  uq  =  0.3  and  uq  =  0.45  (giving  u)p  =  0.38),  and  for  a  =  ±u;j.  Notice 
that  only  the  case  for  gain  at  the  higher  frequency  produces  negative  x'  (x'  ~  —  2,  close  to  the  perfect 
lens  condition  for  a  vacuum  background)  and  x"  =  0  at  a  single  frequency  (uj  =  up).  There  is  thus  no 
fundamental  reason  why  a  metamaterial  with  gain  cannot  achieve  these  conditions  and  hence  it  is  in 
concept  possible  to  achieve  the  perfect  lens  condition  for  an  incident  field  with  frequency  lop  [5]  and  to 
avoid  the  deleterious  affect  of  even  small  loss  [1,  45]. 

The  physical  lens  arrangement  will  dictate  the  image  resolution.  One  suggestion  has  been  an  in¬ 
finitely  wide  slab  having  n  =  —1  [5].  Any  implementation  must  involve  a  finite  width  slab.  With  a 
gain  mechanism,  the  reflectivity  of  the  truncating  boundaries  will  regulate  the  spontaneous  emission 
spectrum.  With  a  thin,  wide  lens,  the  influence  of  the  finite  width  can  be  many  orders  of  magnitude 
smaller  than  the  impact  of  loss  (and  material  mismatch)  on  the  evanescent  field  components  of  the  image 
[45].  Subject  to  the  difficult  task  of  achieving  adequate  gain,  this  paves  the  way  to  far-subwavelength 
imaging.  Perfect  imaging  can  be  achieved  with  a  finite  width  slab  by  placing  the  slab  and  the  object 
within  a  pair  of  coplanar  conductors,  assumed  perfect  at  this  stage.  The  boundary  conditions  on  these 
surfaces,  or  image  theory,  produces  a  periodic  field  solution  in  the  transverse  direction.  With  a  point 
source,  the  periodic  extension  involves  an  infinite  array  of  point  sources.  With  n  =  —  1,  this  infinite 
array  is  perfectly  replicated  at  the  image  plane,  with  the  appropriate  slab  thickness.  Therefore,  in  this 
waveguide  configuration,  with  the  object  field  confined  within  the  transverse  spatial  support  of  the  two 
conductors,  the  object  is  perfectly  replicated  at  the  image  plane.  The  efficacy  of  this  approach  depends 
on  the  influence  of  the  finite  conductivity  of  the  conductor  surfaces.  Good  conductors  are  needed  only 
at  the  wavelength  of  operation,  and  control  of  the  reflectivity  at  the  pump  wavelength  may  be  beneficial 
in  transverse  pumping. 
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Figure  16:  Susceptibilities  \  =  x'  +  x"  from  (15),  with  u>i  =  0.3  and  U2  =  0.45,  for:  (a)  a  =  ojf,  giving  a  loss  resonance  at  uq  and 
gain  at  lo 2,  and  (b)  a  =  —  producing  gain  at  and  loss  at  UJ2  ■ 


With  loss  and  gain  resonances,  and  the  gain  being  at  a  higher  frequency,  a  causal  dispersion  model 
suggests  the  perfect  lens  condition  can  be  approached  in  the  limit  of  ideal  parameters:  n  =  —  1  and  a 
suitable  or  suitably  large  transverse  truncation.  How  effectively  these  conditions  can  be  approached  is 
a  technology  issue  and,  most  fundamentally,  this  hinges  on  whether  an  adequate  gain  mechanism  can 
be  realized.  Any  artificial  material  must  also  have  a  satisfactorily  matched  real  part  of  the  refractive 
index  and  an  adequate  dipole  density. 

4  References 

1.  K.  J.  Webb,  M.  Yang,  D.  W.  Ward,  and  K.  A.  Nelson,  “Metrics  for  negative-refractive-index  materials,”  Phys.  Rev. 
E  70,  035602  (2004). 

2.  M.  Yang  and  K.  J.  Webb,  “Poynting  vector  analysis  of  a  superlens,”  Opt.  Lett.  30,  2382-2384  (2005). 

3.  K.  J.  Webb  and  M.  Yang,  “Subwavelength  imaging  with  a  multilayer  silver  film  structure,”  Opt.  Lett.  31,  2130-2132 

(2006). 

4.  V.  G.  Veselago,  “The  electrodynamics  of  substances  with  simultaneously  negative  values  of  t  and  /r,”  Sov.  Phys. 
Uspekhi  10,  509-514  (1968). 

5.  J.  B.  Pendry,  “Negative  refractive  index  makes  a  perfect  lens,”  Phys.  Rev.  Lett.  85,  39663969  (2000). 

6.  K.  J.  Webb  and  M.  Yang,  “Subwavelength  imaging  with  a  multilayer  silver  film  structure,”  Opt.  Lett.  31,  2130-2132 

(2006). 

7.  P.  A.  Belov,  “Backward  waves  and  negative  refraction  in  uniaxial  dielectrics  with  negative  dielectric  permittivity 
along  the  anisotropy  axis,”  Microwave  Opt.  Technol.  Lett.  37,  259-263  (2003). 

8.  K.  J.  Webb  and  M.  Yang,  “Subwavelength  imaging  with  a  multilayer  silver  film  structure,”  Opt.  Lett.  31,  2130-2132 
(2006). 

9.  S.  A.  Ramakrishna  and  J.  B.  Pendry,  “Removal  of  absorption  and  increase  in  resolution  in  a  near-field  lens  via  optical 
gain,”  Phys.  Rev.  B  67,  201101  (2003). 

10.  B.  Wood,  J.  B.  Pendry,  and  D.  P.  Tsai,  “Directed  subwavelength  imaging  using  a  layered  metal-dielectric  system,” 
Phys.  Rev.  B  74,  115116  (2006). 

11.  P.  A.  Belov  and  C.  R.  Simovski,  “Canalization  of  subwavelength  images  by  electromagnetic  crystals,”  Phys.  Rev.  B 
71,  193105-193108  (2005). 

12.  P.  A.  Belov  and  Y.  Hao,  “Subwavelength  imaging  at  optical  frequencies  using  a  transmission  device  formed  by  a 
periodic  layered  metal-dielectric  structure  operating  in  the  canalization  regime,”  Phys.  Rev.  B  73,  113110-113113 
(2006). 

13.  P.  A.  Belov  and  M.  G.  Silveirinha,  “Resolution  of  subwavelength  transmission  devices  formed  by  a  wire  medium,” 
Phys.  Rev.  E  73,  056607-056615  (2006). 

14.  N.  Engheta,  A.  Salandrino,  and  A.  Alu,  “Circuit  elements  at  optical  frequencies:  nanoinductors,  nanocapacitors,  and 
nanoresistors,”  Phys.  Rev.  Lett.  95,  095504  (2005). 

15.  N.  Engheta,  “Circuits  with  light  at  nanoscales:  optical  nanocircuits  inspired  by  metamaterials,”  Science  317,  1698- 
1702  (2007). 

16.  E.  Shamonina,  V.  Kalinin,  K.  Ringhofer,  and  L.  Solymar,  “Imaging,  compression  and  Poynting  vector  streamlines  for 
negative  permittivity  materials,”  Electron.  Lett.  37,  1243-1244  (2001). 

17.  S.  A.  Ramakrishna,  J.  B.  Pendry,  M.  C.  K.  Wiltshire,  and  W.  J.  Stewart,  “Imaging  the  near  field,”  J.  Mod.  Opt.  50, 
1419-1430  (2003). 

18.  R.  Wangberg,  J.  Elser,  E.  E.  Narimanov,  and  V.  A.  Podolskiy,  “Nonmagnetic  nanocomposites  for  optical  and  infrared 
negative-refractive- index  media,”  J.  Opt.  Soc.  Am.  B  23,  498-505  (2006). 


17 


19.  D.  M.  Pozar,  in  Mircowave  Engineering ,  3rd  ed.  (John  Wiely  &  Sons,  Inc.,  2005),  pp.  412-416. 

20.  Y.-F.  Huang  et  al. ,  “Improved  broadband  and  quasi-omnidirectional  anti-reflection  properties  with  biomimetic  silicon 
nanostructures,”  Nature  Nanotech.  2,  770-774  (2007). 

21.  Handbook  of  Optical  Constants  of  Solids,  E.  D.  Palik,  ed.,  (Academic  Press.,  New  York,  1998). 

22.  K.  J.  Webb  and  J.  Li,  “Analysis  of  transmission  through  small  apertures  in  conducting  films,”  Phys.  Rev.  B  73,  33401 
(2006). 

23.  K.  J.  Webb  and  J.  Li,  “Resonant  Slot  Optical  Guiding  in  Metallic  Nanoparticle  Chains,”  Phys.  Rev.  B  72,  201402R 
(2005). 

24.  K.  J.  Webb  and  J.  Li,  “Resonant  waveguide  field  enhancement  in  dimers,”  Opt.  Lett.  31,  3348-3350  (2006). 

25.  C.  Delerue,  G.  Allan,  and  Y.  M.  Niquet,  “Collective  excitations  in  charged  nanocrystals  and  in  close-packed  arrays 
of  charged  nanocrystals,”  Phys.  Rev.  B  72,  195316  (2005). 

26.  L.  E.  Brus,  J.  Chem.  Phys.  80,  4403  (1984). 

27.  S.  Schmitt-Rink,  D.  A.  B.  Miller,  and  D.  S.  Chemla,  “Theory  of  the  linear  and  nonlinear  optical  properties  of 
semiconductor  microcrystallites,”  Phys.  Rev.  B  35,  8113-8125  (1987). 

28.  P.  G.  Bolcatto  and  C.  R.  Proetto,  “Shape  and  dielectric  mismatch  effects  in  semiconductor  quantum  dots,”  Phys. 
Rev.  B  59,  12487-12498  (1999). 

29.  A.  Franceschetti,  H.  Fu,  L.  W.  Wang,  and  A.  Zunger,  “Many-body  psuedopotential  theory  of  excitons  in  InP  and 
CdSe  quantum  dots,”  Phys.  Rev.  B  60,  1819-1829  (1999). 

30.  Z.  Yu,  J.  Li,  D.  B.  O’Connor,  L.  Wang,  and  P.  F.  Barbara,  “Large  resonant  Stokes  shift  in  CdS  nanocrystals,”  J. 
Phys.  Chem.  B  107,  5670-5674  (2003). 

31.  D.  .J.  Norris,  A.  L.  Efros,  M.  Rosen,  and  M.  G.  Bawendi,  “Size  dependence  of  exciton  fine  structure  in  CdSe  quantum 
dots,”  Phys.  Rev.  B  53,”  16347-16354  (1996). 

32.  M.  Nirmal,  B.  O.  Dabbousi,  M.  G.  Bawendi,  J.  .].  Macklin,  J.  K.  Trautman,  T.  D.  Harris,  and  L.  E.  Brus,  “Fluorescence 
intermittency  in  single  cadmium  selenide  nanocrystals,”  Nature  383,  802-804  (1996). 

33.  C.  B.  Murray,  C.  R.  Kagan,  and  M.  G.  Bawendi,  “Synthesis  and  characterization  of  monodisperse  nanocrystals  and 
close-packed  nanocrystal  assemblies,”  Annu.  Rev.  Mater.  Sci.  30,  545-610  (2000). 

34.  X.  Peng,  .J.  Wickham,  and  A.  P.  Alivisatos,  “Kinetics  of  II- VI  and  III-V  colloidal  semiconductor  nanocrystal  growth: 
focussing  of  size  distributions,”  J.  Am.  Chem.  Soc.  120,  5343-5344  (1998). 

35.  J.  C.  Maxwell-Garnett,  “Colours  in  metal  glasses  and  in  metallic  films,”  Philos.  Trans.  R.  Soc.  Lond.  203,  385-420 
(1904). 

36.  J.  C.  Maxwell-Garnett,  “Colours  in  metal  glasses,  in  metallic  films,  and  in  metallic  solutions.  II,”  Philos.  Trans.  R. 
Soc.  Lond.  205,  237-288  (1906). 

37.  R.  W.  Cohen,  G.  D.  Cody,  M.  D.  Coutts,  and  B.  Abeles,  “Optical  properties  of  granular  silver  and  gold  films,”  Phys. 
Rev.  B  8,  3689-3701  (1973). 

38.  J.  I.  Gittleman  and  B.  Abeles,  “Comparison  of  the  effective  medium  and  the  Maxwell-Garnett  predictions  for  the 
dielectric  constants  of  granular  metals,”  Phys.  Rev.  B  15,  3273-3275  (1977). 

39.  D.  R.  Smith,  S.  Schultz,  P.  Markos,  and  C.  M.  Soukoulis,  “Determination  of  effective  permittivity  and  permeability 
of  metamaterials  from  reflection  and  transmission  coefficients,”  Phys.  Rev.  B  65,  195104  (2002). 

40.  S.  L.  Ray,  “Numerical  Disperion  and  Stability  Characteristics  of  Time-Domain  Methods  on  Nonorthogonal  Meshes,” 
IEEE  Trans.  Antennas  Propagat.  41,  233-235  (1993). 

41.  L.  Brillouin,  Wave  propagation  in  periodic  structures  (Dover,  1953). 

42.  V.  Myroshnychenko  and  C.  Brosseau,  “Finite-element  method  for  calculation  of  the  effective  permittivity  of  random 
inhomogeneous  media,”  Phys.  Rev.  E  71,  016701  (2005). 

43.  L.  Chen,  S.  He,  and  L.  Shen,  “Finite-Size  Effects  of  a  Left-Handed  Material  Slab  on  the  Image  Quality,”  Phys.  Rev. 
Lett.  92,  107404  (2004). 

44.  D.  R.  Smith,  D.  Schurig,  M.  Rosenbluth,  S.  Schultz,  S.  A.  Ramakrishna,  and  J.  B.  Pendry,  “Limitations  on  subd¬ 
iffraction  imaging  with  a  negative  refractive  index  slab,”  Appl.  Phys.  Lett.  82,  1056  (2003). 

45.  M.  Yang  and  K.  J.  Webb,  “Poynting  vector  analysis  of  a  superlens,”  Opt.  Lett.  30,  2382-2384  (2005). 

46.  R.  A.  Shelby,  D.  R.  Smith,  and  S.  Schultz,  “Experimental  verification  of  a  negative  index  of  refraction,”  Science  292, 
77-79  (2001). 

47.  M.  Stockman,  “Criterion  for  negative  refraction  with  low  optical  losses  from  a  fundamental  principle  of  causality,” 
Phys.  Rev.  Lett.  98,  177404  (2007). 

48.  J.  D.  Jackson,  Classical  Electrodynamics,  3  ed.  (Wiley,  New  York,  NY,  1999). 

49.  E.  C.  Titchmarsh,  Introduction  to  the  Theory  of  Fourier  Integrals  (Oxford  University  Press,  Oxford,  UK,  1948). 

50.  L.  Brillouin,  Wave  Propagation  and  Group  Velocity  (Academic  Press,  New  York,  NY,  1960). 

51.  A.  Yariv,  Quantum  Electronics,  2  ed.  (Wiley,  New  York,  1975). 


